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ABSTRACT 

We have worked out simple analytical formulae that accurately approximate the relationship between 
the position of the source with respect to the lens center and the amplification of the images, hence the 
lens cross section, for realistic lens profiles. We find that, for essentially the full range of parameters 
either observationally determined or yielded by numerical simulations, the combination of dark matter 
and star distribution can be very well described, for lens radii relevant to strong lensing, by a simple 
power-law whose slope is very weakly dependent on the parameters characterizing the global matter 
surface density profile and close to isothermal in agreement with direct estimates for individual lens 
galaxies. Our simple treatment allows an easy insight into the role of the different ingredients that 
determine the lens cross section and the distribution of gravitational amplifications. They also ease 
the reconstruction of the lens mass distribution from the observed images and, vice- versa, allow a fast 
application of ray-tracing techniques to model the effect of lensing on a variety of source structures. 
The maximum amplification depends primarily on the source size. Amplifications larger than rj 20 are 
indicative of compact source sizes at high-z, in agreement with expectations if galaxies formed most of 
their stars during the dissipative collapse of cold gas. Our formalism has allowed us to reproduce the 
counts of strongly lensed galaxies found in the H- ATLAS SDP field. While our analysis is focussed 
on spherical lenses, we also discuss the effect of ellipticity and the case of late-type lenses (showing 
why they are much less common, even though late-type galaxies are more numerous). Furthermore 
we discuss the effect of a cluster halo surrounding the early-type lens and of a supermassive black hole 
at its center. 

Subject headings: galaxies: elliptical - galaxies: high redshift - gravitational lensing - submillimeter 

1. INTRODUCTION 

The discovery rate of strong galaxy-scale lens systems has increased dramatically in recent years mostly thanks to 
spectroscopic lens searches and, most recently, to surveys of sub-millimeter galaxies. Spectroscopic searches, such as 
the Sloan Lens Advanced Camera for Surveys (SLACS) survey (Bolton et al. 2006, 2008; Auger et al. 2009) or the 
BOSS (Baryon Oscillation Spectroscopic Survey) Emission-Line Lens Survey (BELLS; Brownstein et al. 2012) or the 
optimal line-of-sight (OLS) lens survey (Willis et al. 2006) or the Sloan WFC (Wide Field Camera) Edge-on Late-type 
Lens Survey (SWELLS; Treu et al. 2011), rely on the detection of multiple background emission lines in the residual 
spectra found after subtracting best-fit galaxy templates to the foreground-galaxy spectrum. 

Sub-millimeter surveys were predicted (Blain 1996; Perrotta et al. 2002, 2003; Negrello et al. 2007), and demon- 
strated (Negrello et al. 2010), to be an especially effective route to efficiently detect strongly lensed galaxies at high 
redshift because the extreme steepness of number counts of unlensed high-z galaxies implies a strong magnification 
bias so that they are easily exceeded by those of strongly lensed galaxies at the bright end. Also, gravitational lensing 
effects are more pronounced for more distant sources. But high-z galaxies are frequently in a dust-enshrouded active 
star formation phase and therefore are more easily detected at far-IR/sub-mm wavelengths, while they are very opti- 
cally faint. Negrello et al. (2007) predicted that about 50% of galaxies with 500 /im flux densities above « 100 mJy 
would be strongly lensed, with the remainder easily identifiable as local galaxies or as radio-loud Active Galactic Nuclei 
(AGNs). This prediction was supported by the mm-wave South Pole Telescope (SPT) counts (Vieira et al. 2010). 
But a spectacular confirmation came with the results of the Herschel Astrophysical Terahertz Large Area SurvejQ 
(H- ATLAS; Eales et al. 2010) for the Science Demonstration Phase (SDP) field covering about 14.4 deg 2 . Five out of 
the 10 extragalactic sources with S^oo M m ^ 100 mJy were found to be strongly lensed high-z galaxies, four are z < 0.1 
spiral galaxies and one is a flat-spectrum radio quasar (Negrello et al. 2010). Gonzalez-Nuevo et al. (2012) presented a 
simple method, the Herschel- ATLAS Lensed Objects Selection (HALOS), aimed at identifying fainter strongly lensed 
galaxies. This method gives the prospect of reaching a surface density of ~ 2 deg -2 for strongly lensed candidates, 
i.e., the detection of ~ 1000 high-z strongly lensed galaxies over the full H- ATLAS survey area (w 550 deg 2 ; Eales et 
al. 2010). 

It should be noted that the (sub-)mm selected lensed galaxies are very faint in the optical, while most foreground 
lenses are passive ellipticals (Auger et al. 2009; Negrello et al. 2010), essentially invisible at sub-mm wavelengths. 
This means that the foreground lens is 'transparent' at (sub-)mm wavelengths, i.e. does not confuse the images of the 
background source. Therefore, the (sub-)mm selection shares with spectroscopic searches the capability of detecting 
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lensing events with small impact parameters and has the advantage that, in most cases, there is no need to subtract 
the lens contribution to recover the source images within the effective radii of the lenses. Also, compared to the optical 
selection, the (sub-)mm selection allows us to probe earlier phases of galaxy evolution which have typically higher 
lensing optical depths. This makes this technique ideal for tracing the mass density profiles of elliptical galaxies over 
a broad redshift range and for probing their evolution with cosmic time. 

Samples of strongly lensed galaxies are further enriched by the, to some extent complementary, imaging surveys 
(Cabanac et al. 2007; Faure et al. 2008; Kubo & Dell' Antonio 2008; Ruff et al. 2011) which look for arc-like 
features, and by radio surveys (Browne et al. 2003). All that holds the promise of a fast increase of the number of 
known strongly lensed sources, fostered by the forthcoming large area optical (e.g., Oguri & Marshall 2010) and radio 
(SKA, Square Kilometer Array) surveys (e.g., Koopmans, Browne, & Jackson 2004). A simple, efficient, analytical 
tool applicable to the analysis of large samples of galaxy-scale lenses is therefore warranted. 

In this paper we work out exact and approximate solutions of the lens equation based on a realistic model for the 
mass density profile of the lens (§[5]) and exploit them to reckon the lensing probability as a function of the source 
redshift (§[3]). As an application of these results, following on the study of the high-z luminosity function of galaxies 
measured by the H- ATLAS survey (Lapi et al. 2011), we compute, in §|4j the number counts of strongly lensed sub-mm 
galaxies implied by the physical model of galaxy formation and evolution formulated by Granato et al. (2004) and 
further developed by Lapi et al. (2006) and Mao et al. (2007). The model counts are compared with the observational 
estimates by Negrello et al. (2010) and by Gonzalez-Nuevo et al. (2012). While our study is focussed on early-type 
lenses, assumed to be circularly symmetric, the cases of late-type lenses, of ellipticity, of the presence of a super-massive 
black hole in the galactic nucleus, and of super-galactic structures are discussed in §[5j Finally, our main results are 
summarized in §[6j 

Throughout the paper we adopt the standard flat ACDM cosmology (see Komatsu et al. 2011) with current matter 
density parameter = 0.27 and Hubble constant H — 72 km s _1 Mpc -1 . 



2. LENSING CROSS SECTION 
2.1. Lens mass models 

We focus here on galaxy-scale lensing, i.e., on those lensing events where the deflector is a single/isolated early- type 
galaxy. The discussion of the effect of the more extended dark matter (DM) halo of a group or cluster in which the 
lens galaxy may reside, or of the disk of later type galaxies is deferred to §[5] 

More specifically, we assume that the lens galaxy is associated to a DM halo of mass Mh in the range 10 n - 4 -10 13 - 5 M Q 
virialized at redshift z^ v > 1.5. The redshift and the lower mass limit are crudely meant to single out galactic halos 
associated with individual spheroidal galaxies. Disk-dominated (and irregular) galaxies are primarily associated with 
halos virializing at z^ v < 1.5, which may have incorporated halos less massive than 10 11 ' 4 Mq virialized at earlier times, 
that form their bulges. The upper mass limit to individual galaxy halos comes from weak- lensing observations (e.g., 
Kochanek & White 2001; Kleinheinrich et al. 2005), kinematic measurements (e.g., Kronawitter et al. 2000; Gerhard 
et al. 2001), and from a theoretical analysis on the velocity dispersion function of spheroidal galaxies (Cirasuolo et 
al. 2005). The same limit is also suggested by modeling of the spheroids mass function (Granato et al. 2004), of the 
quasar luminosity functions (Lapi et al. 2006), and of the sub-mm galaxy number counts (Lapi et al. 2011). 

For the DM we adopt a standard NFW (Navarro, Frenk, & White 1996) profile (e.g., Lokas & Mamon 2001) 

m = JWh f c c 2 

PH(r) 47ri?|(r/ J R H )(l + C r/i?H) 2 ' {> 

where c is the concentration parameter and f c = [log(l + c) — c/(l + c)] . The halo virial radius i?n is given by (e.g., 
Bryan & Norman 1998; Barkana & Loeb 2001) 



i? H = 62 



M H ^ 1/3 
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in terms of the virialization redshift of the lens zi tV , of the evolved density parameter Clyi(z) = J1m (1 + z) 3 /[Q,m (1 + 
z) 3 + 1 — £Im], and of the overdensity threshold for virialization A v (z) = 18 7r 2 + 82 [f2M(z) — 1] — 39 [Om(^) — l] 2 - 
For example, Zg tV — 2.5 and Mh = 10 13 Mq correspond to a virial size i?n ~ 200 kpc. Since the NFW profile yields a 
logarithmically diverging mass, we set the edge of the halo at J?h- 

Hereafter we adopt zg_ v = 2.5 as our fiducial value. This choice is motivated on considering that massive early-type 
galaxies at zt ~ 0.7 feature relatively old ages ^3 — 5 Gyr of their stellar populations, and formed the bulk of their 
stars over a timescale of order 1 Gyr (for a review, see Renzini 2006; also Lapi et al. 2011). These evidences point 
toward a virialization redshift of the host halo in the range z^ jV ~ 1.5 — 3.5, consistent with the distribution of creation 
redshifts found in numerical simulations for the massive halos considered here (e.g., Moreno et al. 2009). Note that for 
an early-type lens the observation redshift (zt ~ 0.7; Gonzalez-Nuevo et al. 2012) is, generally, substantially lower than 
the virialization one (z£ lV ~ 2.5). Therefore the frequently made approximation zg w z^ v leads to a large overestimate 
of the halo size and, indirectly, to an underestimate of the lensing probability. 
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Numerical simulations indicate that the concentration c depends on halo mass and redshift as (Prada et al. 2011) 

(3) 



,. r \ -0.074 



10 13 M C .J \ 1.7 



with a scatter of about 20%. For a lens with Mh = 10 13 M at redshift zg = 0.7 we have c ~ 5, that we adopt as our 
fiducial value. 

The DM to baryon ratio Mh/M* in early-type galaxies is generally in the range 10 — 70. In fact, this quantity can 
be roughly bounded from below by the cosmological DM to baryon mass ratio (see Komatsu et al. 2011) that takes on 
values around 6, and from above by the DM to stellar mass ratio that statistical arguments (see Shankar et al. 2006; 
Lagattuta et al. 2010; Moster et al. 2010) estimate to be around 70. We take Mh/M* = 30 as our fiducial value. 

The impact on the lensing probability of different choices for zg tV , for c [including the use of the halo-mass dependent 
expression of Eq. ([3J], and for Mh/M* is discussed in §[3l 

For the stellar component we adopt the 3-D Sersic profile (Prugniel & Simien 1997) 

P*(r) = -Ml* ( ^V"" e-M'/*^", (4) 



47ri?3 nT(2n) \R e 

where R e is the effective radius, n is the Sersic index, b n = 2n— 1/3 + 0. 009876/n, and a n = 1 — 1.188/2n + 0.22/4n 2 . 
The effective radius R e is related to the stellar mass by (Shen et al. 20033; Hyde & Bernardi 2009) 

M* ^°' 55 

110 



R ° * L28 ( W^) kpc ' (5) 



found to hold (with a w 30% scatter) for Zi<l. Several recent observational studies have found that massive, passively 
evolving galaxies at z > 1 are much more compact than local galaxies of similar stellar mass (Fan et al. 2010 and 
references therein). The study by Maier et al. (2009), with high spectroscopic completeness, finds however that the size 
evolution at fixed mass is modest (~ 25%) from z ~ 0.7 to z = 0, i.e. up to our reference value of the lens redshift. We 
have checked that decreasing R e at fixed M* by a factor f e >l increases somewhat the probability of amplifications 
only in the range 3 < fi < 5; the effect becomes essentially independent of the decrease factor for f e > 1.5. 

Values of the Sersic index for massive early-type galaxies are in the range n tzs 3 — 10, with a tendency for more 
massive systems to feature higher values (e.g., Kormendy et al. 2009). Early-type galaxies with n = 2 generally are 
either dwarf spheroidals or contain a substantial disk component and do not obey Eq. ([3]). We will consider a fiducial 
value n = 4, corresponding to the classical de Vaucouleurs (1948) profile. Again, the effect on the lensing probability 
of different values is investigated in §G3 

To sum up, we adopt a two-component model, made of a stellar component with a Sersic profile plus a DM halo 
with a NFW profile. Given the halo mass Mh and the virialization redshift, zt lV , the total mass distribution of the 
lens galaxy is specified by the parameters zi, M-r/M+, and n [after Eq. ([3]) c is determined by Mh and zg\. Hereafter 
this two-component model will be referred to as the 'SISSA model0. The lensing probability distribution yielded by 
this model will be compared with those yielded by two other commonly used models. The first (hereafter referred to 
as 'NFW model') just consists in adopting a pure NFW profile, hence neglecting the effect of the baryons. The second 
(hereafter referred to as the 'SIS model') adopts the classic singular isothermal sphere density profile 

/ \ Cqr^ 1 Mh / r \ ~ 2 

^ sW = 2S^ = 4^Uj ' (6) 

where ctsis is the 1-D velocity dispersion of the overall mass. The second equality follows from the commonly used 
assumption ctsis ~ Vh/V^, in terms of the halo circular velocity Vjf = G Mtj/Rh. 

2.2. Surface density profile 

The surface density writes 



E(s) = J Ar^==p{r), (7) 

s being the radial coordinate projected on the plane of the sky. It is important to remember that the surface density 
becomes effective for strong lensing when it exceeds the critical threshold 

c 2 (l + z e )D s 



4nG D e D es 

corresponding to a convergence k = 1 for a thin lens. Here D s , Di, and D^ s = D s — D[ are the comoving angular 
distances (also called proper motion distances; see Kochanek 2006) from the source at z s to the observer at z = 0, 

2 Note that there is a typo in the normalization factor given in Table 1 of Shen et al. (2003). The correct value is given in Cimatti et 
al. (2008). 

3 From the name of our main institution. The acronym 'SISSA' and 'SIS' sound as close as the corresponding model results are. 
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from the lens at zt to the observer at z = 0, and from the source at z s to the lens at zt, respectively. In a flat Universe 
the comoving angular distances are defined as = (c/H ) f*' dz [ftjf (1 + z) 3 + 1 — Q.m]^ 1 ^ 2 - The angular diameter 
distance is Dy /(1 + Zj) and the luminosity distance is (1 + zj). 

In Fig. Q]the surface mass density yielded by the SISSA model for a lens at zg = 0.7 with Mh = 10 13 M Q , zg tV — 2.5, 
c = 5, n = 4, and Mh/M* = 30, is compared with those yielded by the Sersic, NFW and SIS laws. The horizontal 
grey line represents the critical surface density when the source is at z s = 2.5 and the lens at zg = 0.7. It shows that, 
in this case, only the matter located at s < 10~ 15 i?H is effective for strong lensing. The stellar contribution to the 
surface density dominates for s < 10 -2 i?n, while the DM contribution takes over at larger radii. 

As illustrated by Fig. [T] in the radial range —2.5 < log(s/i?n) S — 1> which generally contributes most to the 
gravitational deflection, the combination of the stellar and DM components to the total surface density closely mimics 
a power law 

S( S ) = S (A \ (9) 

In Table[TJwe show that, at fixed halo mass, Mh, both the normalization Eo at the reference radius so ~ 10~ 2 i?H, and 
the power-law index 77 0.8 — 0.9 are only weakly dependent on the parameters of the mass distribution. Specifically, 
7/ and £0 slightly increase with decreasing Mh/M* (i.e. for higher stellar contributions in the inner region), with 
increasing halo concentration (i.e. for a higher DM contribution in the inner region), and with decreasing Sersic index 
(i.e. for a stellar distribution more concentrated at the center). At fixed parameters of the lens mass distribution, So 
(obviously) decreases while r\ increases with decreasing halo mass. 

The slopes of the power-law approximation are in good agreement with those determined from the stellar dynamics 
(e.g., Thomas et al. 2011), from the globular cluster/planetary nebulae kinematics (e.g., Rodionov & Athanassoula 
2011), from the HI gas disk/ring (e.g., Weijmans et al. 2008), from the profile of the X-ray emission (e.g., Churazov 
et al. 2010; Humphrey & Buote 2010), and from gravitational lensing in individual galaxies (e.g., Koopmans et al. 
2009; Spiniello et al. 2011; Barnabe et al. 2011). All these observational determinations concur in indicating that the 
overall density profile is roughly isothermal (at least in the inner region), with surface density slopes around 77 « 1. 

In comparing the observational determinations of the slopes of the density profiles it must be kept in mind that the 
profiles are not real power laws. This implies that the slope of the volume density, T}v, is not simply related to the 
slope of the surface density, 77, by rjy = 77 + 1, as in the power-law case. For the combination of Sersic and NFW 
profiles considered here we find that the best fit slopes over the radial range relevant for strong lensing are related 
by ?7y w 77 + 1.2, due to projection effects. Thus, the slightly under-isothermal values of 77 yielded by the SISSA 
model (see Tabled]) are fully consistent with the slightly over-isothermal values of 77^ (in the range 1.9 — 2.3) found 
by Koopmans et al. (2009) and Barnabe et al. (2011). The latter authors also find hints of a steeper slope 77 for the 
least massive systems, consistent with the trend coming out of Table [TJ Finally, note that during the galaxy lifetime 
the mass density profile may evolve under the influence of various physical processes; we address the issue in 



2.3. Lens equation: exact and approximated solutions 

In the following we consider circular lenses, deferring to §[5] the discussion on the effect of ellipticity. In such a case, 
the light rays coming from a distant, point- like source are deflected by an angle 



2 / 9 ,.^( s )_ 4 GM(<0 D 



am,z s )=-- J o M^ = ^^^, (10) 

where 9 is the angular distance of the image from the lens center, $ = s (1 + zi)/Di, and M (< £) = 2ir f~ ds s £(s) is 
the mass within the projected radius £ = 9Di/(l + zt). 

The relation between the position of the source and of its (possibly multiple) lensed image(s) relative to the observer 
is described by the lens equation 

p = e-ja<\e\)- (11) 

Here (3 and 8 are the angles formed by the source and by its images with the optical axis, i.e. with the imaginary line 
connecting the observer and the center of the lens mass distribution. Solving the lens equation means finding all the 
positions 8 of the images for a given source position f3. 
The amplification of the images can be computed as: 

1 da a 

u = with A + = 1 and A_ = 1 . (12) 

A + A_ d8 8 

If either of the two quantities A± vanishes, the magnification formally diverges. Thus the condition X± = defines 
critical curves in the lens plane and corresponding caustics in the source plane. The magnification can be positive 
or negative, implying that the image has positive or negative parity or, equivalently, is direct or reversed. The total 
magnification /x to t is the sum of the absolute values of the magnifications for all the images, i.e., ignoring parity. 

A point-like source perfectly aligned with the observer and the center of the foreground mass distribution is lensed 
into a ring of radius 8e, called Einstein ring. For a circular lens it coincides with the critical curve defined by A_ = 0. 



EFFECTIVE MODELS FOR GRAVITATIONAL LENSING 



5 



The other critical curve defined by A + = 0, if present, approximately coincides with an inner ring of radius Bi < Be- 
In Fig. [2] we show an example of the numerical solutions of the lens equation for the SISSA model. The source is 
located at redshift z s = 2.5 and the lens at redshift zt = 0.7. The parameters of the lens mass distribution are as 
in Fig. [TJ The Einstein ring is located at Be ^ 1.82". In addition, the SISSA model (and also the NFW model, but 
not the SIS model) features another, inner critical curve at 9j ps 0.32". In the top panel we illustrate the relation 
between the angular positions of the source j3 and of the images 9, normalized to the position 9e of the Einstein ring. 
In the bottom panel, we illustrate the relation between the angular position of the source P and the amplification of 
the different images, including the total magnification summed over them. 

As it can be seen from the bottom panel of Fig. 2, the SISSA model features one image for j3 > Pi and three images 
for j3 < Pi, where Pi ps 0.5 9e is the location of the outer caustic corresponding to the inner critical curve at 9i. Of 
the three images, however, the one closest to the optical axis (dot-dashed line) is strongly demagnified, while the other 
two (solid and dashed lines) are amplified by different amounts and feature opposite parity. At the location of the 
outer caustic P = Pi, the second and third images are degenerate with infinite magnification (in the plot the red solid 
curve referring to the total magnification should go to /i —¥ oo at P = Pi, but the divergence has been smoothed out 
for clarity). 

In Fig. [3] we illustrate how the angular positions of the Einstein ring 9e and of the inner critical curve 9i depend on 
the lens mass and redshift. In the top panel we focus on the SISSA model, and show how the position of the Einstein 
ring and of the inner critical curve vary as a function of the lens redshift z&, for different source redshifts z s . The 
bottom panel elucidates the dependence of 9e and 9i on the lens mass for different lens models, at fixed z s — 2.5 and 
zg = 0.7. Plainly, 9e and 9i increase with the halo mass since, after Eq. (fTDj) . more mass implies wider bend angles. 
The SISSA model yields larger values of 9e than the NFW and SIS models but smaller values of 9i than the NFW 
model (the SIS model has no inner critical curve). 

These behaviors, and the comparisons among different mass models, may be better understood looking at approxi- 
mate solutions of the lens equations obtained using the power-law description of the projected mass density [Eq. ©] 
that, as shown by Fig. [TJ approximates quite well a realistic density profile in the radial range most relevant for strong 
lensing. For the SISSA model, r\ ps 0.8 — 0.9 is weakly dependent on the parameters of the mass distribution (see 
Table [1]). The SIS model has a slightly steeper slope (7731s = 1), while, in the range most relevant for strong lensing, 
the NFW can be approximated with a flatter slope ttnfw ~ 0.3 — 0.4. 

Under the power law approximation the deflection angle due to the lens potential within a circle of angular radius 

e is 

a(9) = (13) 

Here and in the following the over-bar (e.g., 9 = 9/9e) denotes normalization to the Einstein radius, that for power-law 
lens models can be simply written as 

l/r, 



>E — CO 



So 



(14) 



with 9 = s (l+zg )/ Dg [cf. Eq. (9)]. TableQ]shows that, at fixed lens and source redshift, the factor [2 E /(2 — if) Ec] 1 /'' 

1 /3 2/3 1 /3 p 

scales approximately as M H ' , implying 9e oc M h ' since trivially 9q cx i?H oc M H ' . This illustrates the power of 
lensing for weighing the halos. Since 

A_ = 1-|0P \+ = l-(l-v)\§\-* (15) 

the magnification of an image is 

^ a7aT = [l-\0\-«] [1-(1- V )\9\-^ (16) 

This equation highlights that, in addition to the critical curve corresponding to the Einstein ring {9 = 9e), for r\ < 1 
there is also an inner ring at 

61 = e E (1 - v) 1/ri - (17) 



To get the magnification from Eq. (|16[) we need to compute the positions, 9, of the images as a function of the angular 
distance P (in units of 9e) of the source from the optical axis by solving the lens equation 

p = 9-l.m-\ (is) 

Summing up over all the images yields the global magnification of the lens model as a function of p. 

For the SIS model (?/ = 1) this gives a constant deflection a = Be, and no radial critical curve. For p > 1, i.e., 
outside Be, the lens equation yields only one image at 9 = 1 + /3 with magnification /i = 1 + 1//3. On the other hand, 
for p < 1, i.e., inside Be, it gives the two images 9± = P ± 1, and related magnifications fi± — 1 ± 1//3; thus the total 
magnification amounts to /1 = |^ + | + = /i+ — /j_ = 2//3. 
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For a generic r\ < 1 it is not possible to solve the lens equation analytically but we find that the numerical solutions 
can be well approximated, over the amplification range 1 < /i < 30 and over the range of parameters explored in 

rr ui m u P =!/(/"- !)'' for /3>/3 7 

Table [TJ by the expressions r_ / \r / _~_ (IQ) 

0= (2/m)" for 

that recover the SIS solutions for 77 = 1. The value fii = r](l — rj)^ 1+1 ^ v corresponds to the location of the inner 
critical curve 9j. 

The goodness of these approximations may be appreciated by eye on looking at the dotted lines in Fig. [2] (both 
panels). They can be useful not only for fast computations of the lensing cross section, but also for constructing 
simulated lensed images of extended sources via the ray-tracing technique or, inversely, for reconstructing the lens 
mass profile from images of lensed sources. 



2.4. Cross sections 

Given the relation between the source position /3 and the total magnification of the images /z, the cumulative cross 
section for lensing magnification, as a function of the lens halo mass Mh and of the lens and source redshifts ze and 
z s , simply writes: 

a(> fj,,M H ,z s ,z e ) = tt/? 2 (^). (20) 

In Fig. @] we illustrate the cumulative lens cross section for the SISSA, NFW and SIS models as a function of Mh, for 
z s = 2.5, Z£ = 0.7 and two values of /x. For massive halos moderate amplifications (/x < 5) are mainly contributed 
by the outer regions of the mass distribution where the DM dominates. Thus the SISSA and the NFW models 
yield very similar values of the cumulative cross section for 2 < fi < 5. On the other hand, strong amplifications 
(fi > 10) are mainly contributed by inner regions where the stellar component becomes important or even dominant. 
Correspondingly, the SISSA cross section is considerably higher than the NFW and close to the SIS one. The different 
shapes of the cross sections as a function of halo mass mainly reflect the different dependencies of 9e on Mh for the 
three models. The value of the projected radius s c at which the surface density equals the critical threshold for lensing, 
S c [Eq. ©], decreases with decreasing Mh, causing a decrease of the lensing cross section. In the case of the NFW 
model, below log(MH/M Q ) = 12 the radius s c lies in the region where the profile is much flatter than that of the 
SIS and SISSA models. As a consequence, it decreases faster with decreasing Mh and, correspondingly, the lens cross 
section drops. 



2.5. Lensing by extended sources 

The above results apply to the idealized case of a point-like source. What about extended sources? The problem 
can be solved via the ray-tracing technique (e.g., Schneider, Ehlers & Falco 1992, p. 304), i.e. by applying the (3 — 6 
and f3 — fi relations to every point of the unlensed light distribution of the source. Some examples are shown in Figs. [5] 
to [7] where z s = 2.5 and ze = 0.7. In the cases of Figs. [5] and [5] the light distribution of the source is modeled as a 
standard Sersic profile with n — 4, while Fig. [7] pertains to a galaxy comprising multiple (4 in this example) luminous, 
giant clumps that appear to be a ubiquitous feature of high-redshift star- forming galaxies (e.g., Forster Schreiber et 
al. 2011). The results are shown for different mass models. 

Figure [5] refers to a zero impact parameter, i.e. the center of the source is aligned with the optical axis. The central 
regions of the source are strongly amplified and deformed into the Einstein ring. The NFW and SISSA models (but 
not the SIS model) yield also an inner critical curve, difficult to see because of the limited resolution of the figure. In 
Fig. [6] the impact parameter is small but non-null. The Einstein ring and the inner critical curve are still present since 
there are regions of the unlensed extended source that are crossed by the optical axis. But the central, brightest region 
is lensed into two deformed and amplified images, one inside and the other outside the Einstein ring. In Fig. [7] the 4 
clumps have small impact parameters (specified in the caption). The lensed image looks like an Einstein ring with a 
knotty structure. High resolution imaging would be necessary to distinguish this configuration from that produced by 
a smooth source profile with an impact parameter close to zero. 

This analysis also provides an useful way to estimate the maximum amplification /i max to be expected from extended 
sources, when their multiple images in the lens plane are not resolved (e.g., Peacock 1982). We compute the quantity 
Mmax as the ratio between the total flux of the lensed images to that of the unlensed source over the image plane (see 
Perrotta et al. 2002 for details). In Fig. [8] we plot /x max versus the offset between the center of the extended source 
and the optical axis for the SISSA (red line), NFW (blue), and SIS (green) lens models. The outcomes are quite 
similar, with the SISSA model providing slightly higher maximum amplifications. For the spherical lenses considered 
here /z max increases monotonically with decreasing offset. 

The figure also illustrates, for the SISSA model, the dependence of /z max on the half stellar mass radius, R e , of the 
source. The dissipative collapse of baryons within the DM halos can result in R e ~ 1 — 3 kpc for log(MH/M©) = 12 — 13 
(Fan et al. 2010). The corresponding values of /i max for close alignments between the source and the lens are in the 
range 30 — 50; /i max depends inversely on the source size, and decreases to ~ 10 for R e = 10 kpc. Thus the magnification 
distribution of sub-mm galaxies can provide information on the scale of the stellar mass distribution for dusty high-z 
galaxies, difficult to determine by other means. We also find that /Lt max increases linearly with the source angular 
diameter distance, which is however very weakly dependent on redshift in the range of interest here (1.5 < z s < 4). 
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3. LENSING PROBABILITIES AND DISTRIBUTIONS 



We compute the lensing optical depth as 
t(z s \ > n) = I 



dzi 



I 



dM H 



d 2 N d 2 V 



<r(> n,Mn,z s ,z e ), 
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where d 2 V/dzd£l is the comoving volume per unit z— interval and solid angle, while d 2 N /dM^dV is the galaxy halo 
mass function (see Shankar et al. 2006 for details) , i.e. the statistics of halos containing one single galaxy. As in Lapi 
et al. (2006) the galaxy halo mass function is computed from the standard Sheth & Tormen (1999, 2002) halo mass 
function by (i) accounting for the possibility that a DM halo contains multiple subhalos each hosting a galaxy and 
(ii) removing halos corresponding to galaxy systems rather than to individual galaxies. We deal with (ii) by simply 
cutting off the halo mass function at a mass of 10 135 M beyond which the probability of having multiple galaxies 
within a halo quickly becomes very high (e.g., Magliocchetti & Porciani 2003). As for (i), we add the subhalo mass 
function, following the procedure described by Vale & Ostriker (2004, 2006) and Shankar et al. (2006), and using the 
fit to the subhalo mass function at various redshifts provided by van den Bosch, Tormen & Giocoli (2005). However, 
we have checked that for the masses and redshifts relevant here (z^ v > 1.5 and 11.4 < log(MH/M ) < 13.5), the total 
(halo + subhalo) mass function differs from the halo mass function by less than 5%. Note that in the computation 
of Eq. (|21[) we do not include the contribution to lensing by massive groups and clusters, for which the parameters of 
the lens mass distribution are different from those adopted here (see discussion in §[5]). 

As illustrated by Fig. [9j the lensing optical depth increases very rapidly with increasing source redshift up to z s ~ 1 , 
grows by a factor of s» 4 between z s — 1 and z s — 2 and by a further factor ps 2.5 at z s — 5. When the magnification 
threshold increases from fi = 2 to \i = 10, t(z s \ > fi) decreases by factors of « 25 — 30 for the SIS and SISSA models 
and by a much larger factor (> 100) for the NFW model. 

The inner integral in Eq. (|2ip gives the lens redshift distribution dp(z£| > fi, z s )/dz£, i.e. the surface density per unit 
redshift interval of lenses located at zt that can produce a strong lensing event with total magnification > (i on a source 
at redshift z s . Examples of lens redshift distributions for the NFW, SIS and SISSA models, with an amplification 
threshold /j, = 2 are shown in Fig. [lOj They are similar with broad peaks at zg 0.6 — 0.7 for z s = 2.5. However 
the SIS and SISSA models yield higher high- z tails than the NFW model. As illustrated for the SISSA model, as the 
source redshift, z s , increases so does the peak of the z^ distribution, the peak broadens and the high-z tail becomes 
more prominent. 

The differentiation of Eq. (|2"Tj) with respect to \x yields (minus) the differential magnification distribution dp(fj,\z s ) /dfi, 
illustrated in Fig. [TT] for the NFW, SIS, and SISSA models with z s — 2.5. The distributions are similar below fi s» 3. 
For higher magnifications the SIS and SISSA models are close to each other and increasingly above the NFW one. 
The only appreciable difference between SIS and SISSA occurs in the range 3 < fi < 6 corresponding to the transition 
between magnifications dominated by the DM to magnifications dominated by the stellar component. As visualized 
in Fig. [11] for the SISSA model, increasing z s increases the normalization of the amplification distribution, reflecting 
the increase in the lensing optical depth, without substantially affecting its shape. 

Figure [T^] elucidates the dependence of the SISSA amplification distribution on the parameters of the lens mass 
model. Lower Mn/M* ratios, i.e. larger amounts of stellar mass, yield higher probability of large amplifications and 
lower values of the transition // between the DM and the star dominated regime (the kink shifts to the left; upper 
left-hand panel). Interestingly, taking into account the mass dependence of the Mji/M* ratio implied by the Lapi et 
al. (2011) galaxy formation model yields results similar to those obtained a constant value Mh/M* = 30. The effect 
of varying the Sersic index of the stellar component (upper right-hand panel) is small for n > 4: as n increases we 
get slightly smaller probabilities of large amplifications and slightly higher probabilities of low amplifications. Note 
that the trend with decreasing n, and especially the feature around \x ~ 5, cannot be extrapolated to lower values of 
n because Eq. ([5]) is no longer appropriate to compute the effective radius R e . In particular n = 1 corresponds to 
exponential disks and disk galaxies have much larger effective radii than given by Eq. ([5]) . Correspondingly, they have 
much smaller surface densities, hence much lower probabilities of large amplifications. The lensing by spiral galaxies 
will be discussed in §[5] 

The bottom left-hand panel shows that higher values of the concentration parameter c of the host DM halo yield 
higher probabilities of large amplifications (coming from the gravitational field in the inner regions of the lens) com- 
pensated by a tiny decrease (imperceptible in the figure) of the probability of small amplifications. Using the mass 
dependent parametrization of Eq. (j3|) gives results almost indistinguishable from those obtained using our fiducial 
value c = 5. Finally, the bottom right-hand panel shows that the probability of large amplifications increases with 
increasing virialization redshift, Zi tV , of the lens, as expected since, for given mass, both the halo radius and the stellar 
effective radius decrease; correspondingly, the surface density increases. As a consequence, adopting zg tV — Zt, as done 
in some analyses, underestimates the probability of large amplifications. 

The distributions of angular separations of the two brightest images for our fiducial lens properties (top panel of 
Fig. [T3|) are very similar for the SIS and SISSA models. In both cases the distribution peaks at 1.58" and has a 
FWHM of 3.58". The NFW model peaks at the same angular separation but has a substantially higher probability 
of smaller angular separations. The dependence on the source redshift is very weak. Note that we are considering 
only galaxy-scale lensing, i.e. we are not taking into account lenses on super-galactic scales. The distribution of the 
differences in the propagation times from the source to the observer for the two brightest images, referred to as time 
delays, yielded by the SISSA model is shifted towards slightly shorter values compared to the SIS model; the NFW 
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model implies still shorter delays, with a significant probability at At < 1 day (for details on the computation of time 
delays see, e.g., Porciani & Madau 2000; Oguri et al. 2002). The dependence on the source redshift is appreciable, 
with the peak of the distribution increasing as a function of z s (see also Oguri et al. 2002; Li & Ostriker 2003). This 
is not in contradiction with the fact that the time delay at given decreases with z s as D s /Di s , since in computing 
the time delay distribution the quantity 1/ Di s must be integrated over as in Eq. (21), and both the outcome of this 
integral and D s increase with z s . 

4. NUMBER COUNTS OF LENSED SUB-MM GALAXIES 

Given the number density d 2 N/dSdz s of unlensed sub-mm galaxies per unit flux density and redshift interval and 
the amplification distribution dp/d[i the observed counts allowing for the effect of lensing can be computed as 



Here we have approximated to unity the factor 1/ < fj, > that would have appeared on the right-hand-side, as appro- 
priate for large-area surveys (see Jain & Lima 2011). For the unlensed counts we adopt the model by Lapi et al. (2011) 
that successfully reproduces the (sub-)mm counts from 250 to 1100 /im. The resulting Euclidean normalized counts 
of point-like sources at 500 /im and 350 /im are shown in Fig. 1141 standard parameters of the lens mass distribution 
have been adopted. The source extension translates into an upper limit on the amplification (Fig. [8]) whose effect on 
the counts is set out in the lower panel of Fig. [15j For compact sources with effective radii in the range R e = 1 — 3 
kpc, expected to be typical for the high-z galaxies of interest here, the predicted counts are quite similar to those for 
point-like sources. The model counts compare quite well with the counts of confirmed strongly lensed galaxies selected 
at 500 /im by Negrello et al. (2010) and with the counts of candidate strongly lensed galaxies selected via the HALOS 
method at 350 /im by Gonzalez-Nuevo et al. (2012). Note that we have corrected the SDP/HALOS counts with the 
flux-dependent purity estimated by the latter authors (see their § 4.3 and Fig. 10), which amounts to about 70% for 
flux densities S'350/jm ^ 85 mJy. 

The completeness of this sample as a function of flux is difficult to assess accurately. The main source of incomplete- 
ness is constituted by the fact that faint (mostly high-redshift) lenses may have been missed by the optical surveys 
exploited by Gonzalez-Nuevo et al. (2012). However, such surveys are expected to be reasonably complete for the 
massive galaxies (M+ >3x 10 11 Mq) acting as lenses in the redshift range zi ~ 0.5 — 0.8 where the lensing probability 
peaks. 

The upper panel of Fig. [15] contrasts SISSA with SIS and NFW models. While even with much larger samples it 
will be hard to discriminate between the SISSA and SIS predictions, the NFW model gives substantially lower counts 
at the brightest flux densities. Thus the counts of bright strongly lensed sources over an H- ATLAS area ten times 
larger than the SDP field (M. Negrello et al., in preparation) can provide a significant test. The lower panel of Fig. [TBI 
illustrates the contributions to the counts of different amplification intervals. Most of the contribution at flux densities 
S350 /im < 80mJy comes from moderate amplifications (2 < fj, < 5). Larger amplifications, up to /i w 30, become 
increasingly important at brighter flux densities. Amplifications /i > 30 are very rare and have little effect on the 
counts, although a few extreme cases may show up in very large area surveys, such as those made by the Planck 
satellite. 

An important application of strong lensing is the possibility of investigating sources fainter than those accessible 
with other means. In particular, sub-millimeter surveys have proven to be an extremely powerful tool to investigate 
the dust enshrouded most active star formation phases of galaxy evolution. However, the detected high-z sub-mm 
sources are generally intrinsically ultra-luminous galaxies, forming stars at extreme rates, (above several hundreds 
solar masses per year; Lapi et al. 2011), and are therefore atypical. Data (Genzel et al. 2006; Foster-Schreiber et al. 
2006) suggest that the most effective star formers in the universe are galaxies with stellar and gas masses ~ 10 11 Mq at 
z w 2 — 3, endowed with star-formation rates of SFR ~ 100 — 200 Mq yr , that generally have sub-mm flux densities 
below the confusion limits of large (sub-)mm surveys (from ground-based single dish telescopes or from Herschcl). 
It is therefore interesting to investigate the flux density range best suited to select such galaxies. Figure \W\ shows 
that, thanks to strong lensing, the H- ATLAS survey will allow us to sample galaxies down to SFR < 200 Mq yr -1 
already at relatively bright flux density limits, where the selection of strongly lensed galaxies is easier (Negrello et 
al. 2010; Gonzalez-Nuevo et al. 2012). For example, we expect sa 0.2 strongly lensed sources per square degree with 
<S350/im ^ 80 mJy and SFR < 200 Mq yr . To find out such galaxies extensive follow-up observations enabling solid 
determinations of the gravitational amplification are necessary. 



The analysis described so far refers to spheroidal lenses. In the case of spiral lenses, there are some significant 
differences. First, the lensing cross-section depends rather strongly on the central surface mass density, which is 
substantially lower in disk compared to spheroidal galaxies. For example, let us consider a spiral galaxy within a halo 
as massive as that of our reference early-type lens (Mh = 10 13 M©), and with the same stellar to DM mass ratio 
30 (this ratio is found to be essentially independent of the Hubble type for massive galaxies; see Mandelbaum et al. 
2006; Shankar et al. 2006). Using the average relation by Tonini et al. (2006b) for the corresponding stellar mass, 




(22) 



5. DISCUSSION 



5.1. Late-type lenses 
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M* ss 3 x 10 11 M , we find a scale radius of the exponential disk profile [E(r) = M ir /{2nR 2 D e~ r / R °)] R D « 7 kpc. 
The interstellar medium adds little: its mass is a small fraction of the stellar mass and it is more diffuse, with a scale 
radius typically of 3 Rd ■ 

In the radial range r < 10~ 15 i?H (with i?n ~ 200 kpc) relevant for strong lensing the surface density of the disk 
is as flat as that of the DM, and is a factor 4 smaller. As a consequence, the lensing properties of late-type galaxies 
are quite close to those of a pure NFW configuration. For comparison, the stellar surface density of an early-type 
galaxy with the same halo mass increases steeply inwards (cf. Fig. [T]); it is a factor sa 3 smaller than the DM's at 
r < 10 _1 ' 5 i?H but a factor 5 larger at r < 10 _2 5 i?H- After Fig. H] this implies that the cross sections (and relatedly 
the number) of massive spiral lens yielding amplifications larger than 10 is lower than that of an early-type with the 
same mass by a factor greater than 3. 

Moreover spheroids, even though less numerous, are, on average, substantially more massive (galaxies with Mh 
as large as 10 13 Mq are much more frequently spheroidals than spirals) and contain the major share of stellar mass 
(Baldry et al. 2004). Furthermore, disk galaxies have generally younger stellar populations than spheroidal galaxies 
(Bernardi et al. 2010; their Fig. 10), indicative of a formation rcdshift < 1; thus they are likely increasingly rarer than 
spheroidal galaxies at substantial rcdshifts. This explains why the contribution of spiral galaxies to the lensed counts 
is subdominant; preliminary evidences of such an expectation come from lens searches with different selection criteria 
(e.g., Auger et al. 2009; Gonzalez-Nuevo et al. 2012). 

5.2. Effect of ellipticity 

Realistic lens models include some ellipticity. In fact, accounting for ellipticity is essential in order to successfully 
reproduce image numbers, image positions and extended lensed images, in most of the observed lensing systems. As 
an example, axially symmetric lenses cannot produce an Einstein cross because the tangential caustic is collapsed to a 
point. Therefore no more than two (in the case of a singularity at the center like in the SIS model) or three (in the case 
of a finite core density like in the NFW or in the SISSA model) images can form. But when ellipticity is non-null, the 
tangential caustic has a finite shape and once the source encompasses it a new pair of images can form (see Fig. [i~7)) . 

A direct comparison between the spherical and the ellipsoidal case can be most easily done considering a Singular 
Isothermal Ellipsoid (SIE hereafter). As we are mostly interested on the lensing statistics we focus on the effect of 
ellipticity on the cross-section for total magnifications cr(/xtot)- We consider the case of a SIE with ellipticity c = 0.2 
and of a SIE with ellipticity c = 0.4, and assume a lens at redshift zg — 0.7, with velocity dispersion a = 350 km s _1 , 
and a source redshift z s = 2.0. The exact choice of these values is irrelevant as one can always work in normalized 
units, using as a reference scale the value of the Einstein radius for a SIS with the same parameters (in this case 
Be ~ 1.85"). We use the public code GLAFld^, a software developed for studying strong gravitational lenses (Oguri 
2010), to solve the lens equation for the SIE model and to estimate the total amplification as a function of the source 
position in the source plan^EL 

On the top panels of of Fig. [17] we show the map of the total amplification as a function of the source position from 
the center of the lens for the SIE model with e = 0.2 (central panel) and for the SIE model with c = 0.4 (right panel). 
Contours of equal amplifications ranging from 2 to 10 are shown in orange. The red curves mark the caustics. For 
comparison, the case of a SIS model with the same lens and source parameters is shown on the top left panel of the 
same figure. 

From these images the cross-section a(fitot) is easily derived and the results are presented in the bottom panel of 
Fig. [17] We see that for amplifications below ~ 3 — 4, the SIE model is almost indistinguishable from the SIS model 
in terms of cross-section. The effect of ellipticity is a 'squeezing' of the equal amplification contours along the major 
axis of the lens (oriented North-South) while leaving the enclosed area almost unaffected. On the other hand, for 
amplifications close to fitot ~ 8 for e = 0.2, and to /Xtot ~ 4 for e = 0.4, the cross-section of the SIE model deviates 
appreciably from the regular circular shape observed for the SIS model, as the source is now close to the inner cross-like 
caustic, and becomes correspondingly smaller (by ~ 30%). The situation is reversed for higher amplifications as the 
cross-section for the SIS model shrinks to a point while that of the SIE model converges to the finite inner caustic. In 
this case the SIE model yields a cross-section that is ~ 50% higher than that given by the SIS model and approaches 
the SIS limit asymptotically for fi to t 3> 100. The exact value of the amplification at which the transition of the 
cross-section from the 'sub-SIS' to the 'super-SIS' regime occurs depends on the adopted value of the ellipticity (it is 
around fitot ~ 20 for e = 0.2 and close to fitot ~ 8 for e = 0.4). In fact, as the ellipticity is increased, the inner caustic 
becomes more extended and the regions of low amplifications in the source plane are consequently more affected. 

In conclusion, compared to the case of a spherical lens, the effect of ellipticity is to slightly decrease (by about 30%) 
the cross-section in a range of amplifications //tot ^4 — 20 (the exact interval depending on the value of the ellipticity) 
and to increase the cross-section (up to 50%) for higher amplifications. 

5.3. Evolution of the mass density profiles 

4 http://www.slac. stanford.edu/~oguri/glafic/ 

5 More precisely we have used the 'mockl' command to randomly populate the source plane with 500000 sources, to solve the lens 
equation and get the amplifications for the individual images for each source position. We have then constructed a map of the total 
amplifications by grouping the source positions into pixels of 0.01" in size. We have finally used such a map to derive the cross section 
cr(/itot). 



10 



A. LAPI ET AL. 



The mass density profile may evolve during the galaxy lifetime under the action of several processes. For example, 
in the early stages of galaxy formation when gas and stars condense toward the inner regions of the system, 'halo 
contraction' may lead to a steepening of the initial, NFW-like DM profile (see Blumenthal et al. 1986; Mo et al. 1998; 
Gnedin et al. 2004). The strength of the effect is widely debated, but recent numerical experiments (see Abadi et al. 
2010; Pedrosa et al. 2010; see also Gnedin et al. 2012) suggest that the classic treatments based on adiabatic invariants 
are likely to be extreme, and that actually in the inner regions the contraction may be inefficient and the density shape 
hardly modified. We have checked that the projected surface density of an overall configuration constituted by a 
Sersic profile for the baryons and a contracted NFW profile for the DM is still well approximated in the radial range 
—2 < log(s/i?n) ;$ — 1 by a power-law shape, with index steeper than the basic SISSA model. Specifically, for a 
lens galaxy at zg — 0.7 with halo mass Mh = 10 13 M Q and fiducial parameters of the mass distribution, we find that 
the slope r\ ps 0.823 of the basic SISSA model is steepened to the value 0.958 according to the prescription for halo 
contraction by Blumenthal et al. (1986), to 0.920 according to that by Gnedin et al. (2004), and to 0.877 according 
to that by Abadi et al. (2010). 

On the other hand, a flattening of the mass density profile may be caused by transfer of energy and/or angular 
momentum from (baryonic and DM) clumps to the DM field during the galaxy formation process (e.g., El-Zant et 
al. 2001; Tonini et al. 2006a). Moreover, at the formation of a spheroid, central starbursts and accretion onto a 
supermassive black hole may easily discharge enough energy with sufficient coupling to blow most of the gaseous 
baryonic mass within the star-forming region out of the inner gravitational well. This will cause an expansion ('puffing 
up') of the stellar and of the DM distributions (see Fan et al. 2010; Ragone-Figueroa et al. 2012), so as to flatten the 
inner profiles. 

A reliable assessment of these competing processes is still lacking, and beyond the scope of the present paper. 
Anyway, it is interesting to point out that, since early- type galaxies generally do not show signs of dissipation at z < 1, 
the steepening of the density profiles due to halo contraction should show up at higher z, though the trend may be 
partly offset by other processes such as the energy transfer or the puffing up mentioned above. 

Recent results by Ruff et al. (2011) and Bolton et al. (2012) show preliminary evidence for a mild evolution in the 
opposite direction, i.e., towards steeper mass profiles at later cosmic times. In fact, the average density slopes steepen 
from values r\ ps 2 at z ps 0.6 toward r\ ps 2.2 at z ps 0.2. The significativity of the detection is still under debate, 
since the trend can be partly explained in terms of variations in the lensing measurement aperture with redshift, 
which favours the sampling of inner, steeper portions of the mass distribution at decreasing z. In any case, we stress 
that around zg ps 0.6 where the lens redshift distribution peaks (see Fig. 10), the measured average density slope 
r)v ~ 2 (corresponding to a surface density slope r\ ps 0.8, see § 2.2) is in excellent agreement with the SISSA model 
outcomes. In addition, even if the trend toward steeper slopes r\ > 2.2 (i.e., ?7y > 1) at z ps 0.2 will be confirmed, our 
lensing analysis based on power-law representations of the surface density would still apply. The overall effect on the 
amplification distribution would be small since the lens redshift distribution is steeply declining for z < 0.5, although 
in the phenomenology of individual lensing systems (e.g., Einstein ring size; presence/absence of inner critical curve) 
the local lenses would tend to behave as ideal SIS more than high-redshift ones. 

5.4. Super-massive black holes in the lens centers 

Super-massive black holes (BHs) are ubiquitous in the nuclei of early-type galaxies (see Ferrarese & Ford 2005 for an 
exhaustive review) . What is the effect of a supermassive BH in the nucleus of a lens galaxy? For an AGN in the source 
the analysis of point source lensing f? !2.4l and 12 . 3[) applies. A thorough discussion will be presented in a subsequent 
paper. 

By itself, an isolated point mass M, features a surface density £(0) oc M, 5d(9)/9 in terms of the Dirac delta function 
6d(0), which yields a deflection profile a(9) oc Mm/9; in our formalism based on Eqs. (9) and (13) this configuration 
corresponds, approximately for £ oc 6^ and exactly for a oc 9 71 , to the limiting powerlaw index r\ — > 2. Then after 
Eq. (18) it is seen that two images are produced independently of the impact parameter (actually there is a third 
image but it has zero magnification), and no inner critical curve is present [cf. § I2.3| . Normalizing angles in units of 

6e — > [Mm (1 + zg) 2 /tt Df Ec] 1 / 2 , the locations of the images are 9 = (/? ± v4 + /3 2 )/2 and their magnifications are 
fi± = ± (/3/v/4 + /3 2 + yjl + 4/j3 2 ± 2)/4 so that the total magnification reads (i = (2 + /3 2 )/^\/4 + P 2 - When the 
source approaches the optical axis at (3 — > 0, the total magnification diverges as /3 _1 and the image positions tend to 
8e (for details, see Kochanek 2006). 

Now we turn to discuss the effect of a supermassive BH at the center of a galactic structure. For sake of definiteness, 
let us consider a lens galaxy with Mh = 10 13 Mq and a mass density profile described by our fiducial SISSA model 
with concentration c = 5, Sersic index n = 4, and dark to baryon mass ratio Mh/M* = 30. To this configuration we 
add a central supermassive BH, with standard BH to stellar mass ratio M./M* ps 2.5 x 10~ 3 . 

The overall surface density is illustrated in Fig. [T8j here, to avoid dealing with the central singularity associated to 
the BH, we rendered its surface density with the powerlaw £(#) oc 9~ n in the limit rj — > 2 (see above). It is seen that 
the BH dominates only at radii s/Rn < 10~ 3 . 

In Fig. [19] we show the solutions of the lensing equation, in terms of the f3 — 9 and /3 — [i relations. These figures 
should be contrasted with Fig. [T] and [5] that represent the same lens configuration without central supermassive BH. 
All in all, the presence of the BH produces three effects: first, the radius of the Einstein ring is increased, although 
only slightly given the small BH contribution to the enclosed mass; second, the inner critical curve is erased, since the 
steep surface density of the point-mass lens dominates the overall behavior for 9 — > 0; third, the central demagnified 
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image can be accompanied by a second detectable central image when the BH mass falls in the range 10 7 5 — 10 8 5 Mq 
(Rusin, Keeton & Winn 2005). 

In Fig. [20] we show that the central supermassive BH has only a minor effect on the differential amplification 
distribution. In this example, we have adopted two values of the BH to stellar mass ratio, and specifically a standard 
one M./M+ — 2.5 x 10~ 3 and a high one 2.5 x 10~ 2 as measured recently in two giant ellipticals by McConnell et 
al. (2011). We have also checked that the mass/redshift dependent M, — M* relationship from the galaxy formation 
model of Lapi et al. (2006, 2011) yields outcomes similar to the former case. 

5.5. Effect of a galaxy cluster halo surrounding the lens galaxy 

The presence of a galaxy cluster halo around an early-type lens can strongly affect the lensing phenomenology of 
individual objects. However, if the lens system is roughly centered on the cluster halo, we find that the effects of 
such events on the amplification distribution are minor. This can be understood on considering the low abundance of 
cluster-sized halos in the halo mass function, rapidly decreasing with increasing redshift, and taking into account that 
the distribution of lens redshifts typically peaks at substantial redshifts (Fig. [T0|) . For example, consider an early-type 
lens with Mr = 10 13 Mq at the center of a cluster halo of 10 15 M©. In the radial range relevant for strong lensing 
(r ~ 10~ 2 - 5 - 10 _1 i? H with a galactic halo size i? H ~ 200 kpc) the cluster halo contributes «2x 10 10 - 7 x 10 9 M 
kpc -2 to the surface density, i.e., an amount comparable or slightly larger than that due to the mass within the galaxy. 
On the other hand, atzx 0.7 an halo of 10 15 Mq is rarer than a halo of 10 13 Mq by a factor 10 6 , implying that the 
contribution of these events to the amplification distribution are negligible. 

On the other hand, some contribution may arise from events in which the lens system is not centered on the cluster 
halo, but is in the vicinity of a galaxy cluster/group in projection. This is because in such instances the cluster/group 
induces an external shear on the lens system, that breaks the spherical symmetry and leads to astroid caustics similar 
to those arising from a non-null ellipticity. As discussed in § 5.2, the lensing cross section is correspondingly affected 
and, considering that the statistics of such events may be higher than that of lens systems centered with the cluster 
halo (though it is difficult to provide an educated estimate), the contribution to the amplification distribution may 
also be non-negligible. 

6. CONCLUSIONS 

In view of the large samples of strongly lensed galaxies that are being/will be provided by large area sub-mm 
(Serjeant 2011; Negrello et al. 2010; Gonzalez-Nuevo et al. 2012), optical (e.g., Oguri & Marshall 2010) and radio 
(SKA) surveys (e.g., Koopmans et al. 2004) we have worked out simple analytical formulae that accurately approximate 
the relationship between the position of the source with respect to the lens center and the amplification of the images 
and hence the cross section for lensing (see Fig. [2]). The approximate relationships are based on a lens matter density 
profile appropriate for early-type galaxies, that comprise most of the lenses found with different selection criteria. 
The adopted profile is a combination of a Sersic profile, describing the distribution of stars, with a NFW profile for 
the dark matter. We find that, for essentially the full range of parameters either observationally determined (for the 
Sersic profile) or yielded by numerical simulations (for the NFW profile) , the combination can be very well described, 
for lens radii relevant for strong lensing, by a simple power law. Remarkably, the power law slope is very weakly 
dependent on the parameters characterizing the matter distribution of the lens (the dark matter to stellar mass ratio, 
the Sersic index, the concentration of NFW profile). For the most common parameter choices, the slope is slightly sub- 
isothermal if we consider the projected profile and slightly super-isothermal if we consider the 3-dimensional profile, 
in good agreement with the results of detailed studies of individual lens galaxies (e.g. Koopmans et al. 2009; Spiniello 
et al. 2011; Barnabe et al. 2011; Ruff et al. 2011; Grillo 2012; Bolton et al. 2012). Our approach implies slightly 
steeper slopes of the total matter density profile for the least massive systems (see Table [T|) ; evidence in this direction 
has been reported by Barnabe et al. (2011). 

Table [1] shows that, if the source and lens redshifts are measured and the halo mass of the lens is reliably estimated, 
the factor [2 E /(2 - 77) £ c ] 1/t? , and hence 9e [see Eq. (|14|)]. varies by no more than 20 — 30% for conceivable variations 
of the parameters of the lens mass distribution. Such small variance paves the way to the possibility of exploiting 
gravitational lensing as a probe of cosmological parameters (Grillo et al. 2008). 

Our simple analytic solutions provide an easy insight into the role of the different ingredients that determine the lens 
cross section and the distribution of gravitational amplifications. The maximum amplification depends primarily on 
the source size. Amplifications larger than w 20, as found for some sub-mm and optical sources (Belokurov et al. 2007; 
Negrello et al. 2010; Swinbank et al. 2010; Brownstein et al. 2012), are indicative of compact source sizes at high-z, 
in agreement with expectations if most of the stars formed during dissipative collapse of cold gas. Similarly, analytic 
formulae highlight in a transparent way the role of parameters characterizing the lens mass profile (Afu, Mr/M* ratio, 
concentration of the DM component, Sersic index of the stellar component), and of the source and lens redshifts. They 
also allow a fast application of ray-tracing techniques to model the effect of lensing on a variety of source structures. 
We have investigated, in particular, the cases of a point-like or of an extended source with a smooth profile, and of a 
source comprising various emitting clumps (as frequently found for high-z active star- forming galaxies). Our formalism 
has allowed us to reproduce the counts of strongly lensed galaxies found in the H- ATLAS SDP field. 

While our analysis is focussed on spherical lenses, we have also discussed the case of disk galaxies (showing why they 
are much less common, even though late-type galaxies are more numerous) and the effect of ellipticity. Furthermore 
we have discussed the effect of a cluster halo surrounding the early-type lens and of a supermassive BH at its center. 
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Figure 1. Surface mass density of an early-type lens. (Blue) dashed line: NFW dark matter profile with Mjj = 10 13 Mq and concentration 
parameter c = 5; (cyan) dot-dashed line: Sersic profile (n = 4) of the stellar component in the proportion Mji/Af* = 30; (red) solid line: 
SISSA model constituted by the sum of the two contributions; (green) triple-dot-dashed line: classical SIS model for the same DM mass; 
(magenta) dotted line: power law relation S(s) oc s~ °- 8 , that provides a good approximation to the SISSA model in the radial range 
relevant for strong gravitational lensing (see text for details). The horizontal grey line represents the critical density for lensing for a source 
redshift z a = 2.5 and a lens redshift zg = 0.7. The projected radius s is normalized to the halo virial radius. 
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Figure 2. Solutions of the lens equation for the SISSA model. Top panel: relation between the angular positions of the source and 
of the images 9, normalized to the position Be of the Einstein ring (the red portion of the curve is zoomed in the inset). Bottom panel: 
relation between the angular position of the source /3 and the amplification /i; of the three images (the black solid, dashed and dot-dashed 
curves refer to different images). The red solid curve is the total magnification. In this example, the source is at z s = 2.5 and the lens at 
%i = 0.7. The parameters of the lens mass distribution are the same as in Fig. [l] The Einstein ring is located at 8g PS 1.82". Note that 
the SISSA and NFW models (but not the SIS model) also feature an inner criti cal c urve at 0i 0.32". In both panels the orange dotted 
curves show the results obtained using the approximate solutions presented in § 12.31 the agreement with numerical solutions is strikingly 
good. 





Figure 3. Top panel: angular positions of the Einstein ring 9e (blue lines) and of the inner critical curve 9j (red lines) for the SISSA 
model, as a function of the lens redshift zg and for different source redshifts: z s = 1.5 (dot-dashed lines); 2.5 (solid lines); and 3.5 (dashed 
lines). The parameters of the lens mass distribution are the same as in Fig. [T] Bottom panel: angular positions of the Einstein ring 9g 
(solid lines) and of the inner critical curve 9j (dashed lines) as a function of the lens mass, for the SISSA (red lines), NFW (blue lines), 
and SIS (green lines) models. The source is at redshift z s = 2.5, the lens at z e = 0.7. 




Figure 4. Lens cross section as a function of the lens mass, for two different amplification thresholds (fi > 2, solid lines, and n > 10, 
dashed lines), and for the SISSA (red lines), NFW (blue lines), and SIS (green lines) models. The source and the lens are at the same 
redshifts as in Fig. \3\ 
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Figure 5. Ray-tracing simulation of the gravitational lcnsing for an extended source at z s = 2.5 with a Sersic profile (n = 4) and effective 
(half light) angular radius 6 S = 0.12", and a lens at zg = 0.7. The lens parameters are the same as in Fig.[T] The impact parameter 8f, (i.e. 
the angular separation between the source and the optical axis) is null in this example. The unlensed source is shown in the top left panel, 
while the lensed image for the SISSA model is in the top right panel, for the NFW model in the bottom left panel, and for the SIS model 
in the bottom right panel. In all panels the origin of coordinates marks the position of the optical axis, while the color scale represents in 
logarithmic units the surface brightness relative to the central value of the unlensed image. Note that the central pixels in the SISSA and 
NFW panels correspond to the strongly demagnified image, and appear in the figure owing to the finite resolution of the simulations. 
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Figure 6. As in the previous figure, but for an impact parameter 8t = 0.35". 
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Figure 7. As in the previous figures, but for four bright spots with impact parameters 0^ = 0.3,0.4,0.15,0.2", respectively. 
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Figure 8. Maximum amplification for an extended source as a function of the offset between the center of the source and the optical axis, 
for the SISSA (red line), NFW (blue), and SIS (green) models. The lens parameters and the source redshift are the same as in Fig.[T] The 
source surface brightness profile is described by a Sersic law with n = 4 and effective radius R e = 3 kpc (solid lines). For the SISSA model 
we also show the results with R e = lkpc (dashed line) and R e = 10 kpc (dotted line). 




Figure 9. Lensing optical depth as a function of the source redshift, for two different amplificationx thresholds (fi > 2, solid lines and 
fi > 10, dashed lines), and for the SISSA (red lines), NFW (blue lines), and SIS (green lines) models. 
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Figure 10. Lens redshift distributions yielded by the SISSA (red lines), NFW (blue line), and SIS (green line) models for a source redshift 
z s = 2.5 (solid lines). For the SISSA model we also show the results for z s = 1.5 (dot-dashed line) and z B = 3.5 (dashed line). 




Figure 11. Amplification distributions yielded by the SISSA (red lines), NFW (blue line), and SIS (green line) models for z a = 2.5 (solid 
lines). For the SISSA model we also show the results for z s = 1.5 (dot-dashed line) and z s = 3.5 (dashed line). 
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Figure 12. Dependence of the amplification distribution for the SISSA model on various parameters: DM to stellar mass ratio Mn/M+ 
(top left panel); Sersic index n (top right panel); DM concentration c (bottom left panel), and lens virialization redshift Z£ v (bottom right 
panel). 
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Figure 13. Top panel: distribution of the separations of the two brightest images of a source at z s = 2.5 for the SISSA (red lines), NFW 
(blue line), and SIS (green line) models. For the SISSA models we also show the results for z s = 1.5 (dot-dashed line) and z s = 3.5 (dashed 
line). Bottom panel: delay time distribution; linestyles are as above. 
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Figure 14. Top panel: Euclidean normalized counts at 500 fira. The open squares and the filled stars represent the total Herschel- 
ATLAS SDP counts (Clements et al. 2010) and the counts of the bright SDP lenses spectroscopically confirmed by Negrello et al. (2010), 
respectively. The solid line illustrates the total model counts comprising the contributions of unlensed (blue dot-dashed line; Lapi et 
al. 2011) and strongly lensed (purple dotted line; SISSA model from the present work) proto-spheroidal galaxies, normal late-type plus 
starburst galaxies (green dashed line; Negrello et al. 2007), radio sources (magenta triple-dot-dashed line; De Zotti et al. 2005). Bottom 
panel: same but at 350 jim. Here the filled squares refer to counts of the SDP candidate lenses selected by Gonzalez-Nuevo et al. (2012). 
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Figure 15. Top panel: Euclidean normalized counts of lensed galaxies at 350 fim. The filled squares refer to the SDP candidate strongly 
lensed galaxies selected by Gonzalez-Nuevo et al. (2012). Our model predictions are also plotted, for the SISSA (red lines), NFW (blue 
lines), and SIS (green lines) models. In the inset the corresponding integral counts are plotted. Bottom panel: same as above, but for 
the SISSA model with different cuts in maximum amplification: /imax = 5, 10, 30, 50, and 100 (green, cyan, blue, magenta, and red line, 
respectively). The line with /i ma i = 100 is indistinguishable from that without amplification cut. 
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Figure 16. Top panel: Euclidean normalized counts of lensed galaxies at 350 /im. Our model predictions on the counts of lensed proto- 
spheroidal galaxies for the SISSA model are plotted for different upper limits to the star formation rate: SFR < 200, 300, 500, and 1000 Af© 
yr - 1 (magenta, blue, green, and red lines, respectively); the latter line is undistinguishable from that without SFR limit. In the inset the 
corresponding integral counts are plotted. Bottom panel: ratio between the counts of lensed galaxies with SFRs below the above thresholds 
and those without any SFR limit. 
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Figure 17. Top panels: map of the total amplification as a function of the source position from the center of the lensing mass, for a SIS 
lens (left), and for a SIE lens with cllipticity e = 0.2 (middle) and e = 0.4 (right). The lens and the source are assumed to be at redshift 
0.7 and 2.5, respectively. The orange lines represent contours of equal amplifications ranging from 2 to 10, while the red point/curve mark 
the caustics. Bottom panel: corresponding cross sections for lensing as a function of the total amplification for the SIS lens (black dashed 
line), and for the SIE lens with e = 0.2 (red solid line) and e = 0.4 (red dot-dashed line). 
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Figure 18. Surface density profile of a configuration constituted by an early-type galaxies and a central super-massive BH. This figure 
should be contrasted with Fig. [l] which represents the same configuration without central super-massive BH. (Blue) dashed line: dark 
matter component with mass Mg = 10 13 Mq and NFW profile with concentration parameter c = 5. (Cyan) dot-dashed line: stellar 
component in the proportion Mjr/M* = 30 relative to the DM with a Sersic n = 4 profile. (Orange) triple-dot-dashed line: super-massive 
BH with mass ratio Af./M* = 2.5 X 1CT 3 . (Red) solid line: SISSA+BH model, sum of the three contributions. Note that here, to avoid 
representing the central singularity, we rendered the point mass surface density with the powerlaw £(#) <x 8~ ri in the limit 77 — >■ 2 (see 
§ 5.4). 




Figure 19. Solutions of the lensing equation for the SISSA model, including a central super-massive BH with mass ratio M./M* = 
2.5 X 10 -3 . This figure should be contrasted with Fig. [2] which corresponds to the same configuration without central super-massive BH. 
Linestyles as in Fig. [2] 
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Figure 20. Differential amplification distribution for a configuration constituted by an early-type galaxies and a central super-massive 
BH. Blue line: SISSA model with fiducial parameter values, without central super-massive BH. Red line: SISSA model plus a super- 
massive BH with standard mass ratio M./M* = 2.5 X 10~ 3 . Green line: SISSA model plus a super-massive BH with a high mass ratio 
M./M* = 2.5 X 1CT 2 . 
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Table 1 

Power law fits of the surface density profile 





Lens parameter 
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Note. — Power law fits of the surface density profile [Eq. {9} over the radial range 

10-2-5 < s //j H < 10-1 for halo masses Mn = 10 , 
10 12 , and 10 13 Mq. The normalization So (in Mq kpc -2 ) refers to a projected radius so = 10 — 2 ijpj- The table shows the dependence 
of the fits on the parameters of the lens mass distribution, namely, dark matter to star mass ratio Mn/M*, Sersic index n, and halo 
concentration c. The case L+ll takes into account the dependence of ratio Mn/M* on Afjj as implied by the Lapi et al. (2011) model. 



